home *** CD-ROM | disk | FTP | other *** search
- #include "yagi.h"
-
- void mprove(double **a,double **alud,int n,int *indx,double *b,double *x)
- {
- int j,i;
- double sdp;
- double *r;
-
- r=dvector(1,n);
- for (i=1;i<=n;i++) {
- sdp = -b[i];
- for (j=1;j<=n;j++) sdp += a[i][j]*x[j];
- r[i]=sdp;
- }
- lubksb(alud,n,indx,r);
- for (i=1;i<=n;i++) x[i] -= r[i];
- free_dvector(r,1,n);
- }
-